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A  COMPARATIVE  STUDY  OF  DIRECTION  FINDING  PERFORMANCE 


by 

Myl^ne  Toulgoat  and  Ross  M.  Turner 


ABSTRACT 


This  report  describes  simulation  studies  of  fast  projection  techniques  for  direction  finding 
of  jammer  signals.  These  techniques  are  based  on  a  data  selection  criterion  applied  directly  to 
the  data  vectors  obtained  at  the  output  of  an  antenna  array.  The  selected  data  vectors  are  used 
in  fast  projection  algorithms  to  estimate  the  jamming  signal  subspace  and  then  to  calculate  the 
directional  spectrum.  Three  such  algorithms  are  described  ranging  from  the  fastest  but  least 
effective  to  the  most  computationally  demanding  but  most  effective.  These  algorithms  provide 
a  better  trade-off  of  performance  versus  computational  load  than  has  heretofore  been  available. 
Relative  performance  is  compared  with  that  of  the  MUSIC  technique  for  test  cases  which 
evaluate  the  minimum  resolution  threshold  achievable  and  the  effect  of  the  jammer  strength 
relative  to  the  receiver  noise. 


RESUME 


Ce  rapport  est  une  etude,  bas^e  sur  des  simulations,  de  la  performance  relative  des 
techniques  de  projection  rapides  pour  trouver  la  direction  de  brouilleurs.  Ces  algorithm.es 
utilisent  un  critfere  pour  la  selection  de  bons  vecteurs  de  donn^es  ^  la  sortie  d’une  antenne 
r^seau.  Les  vecteurs  sdlectionn6s  sont  utilises  dans  des  algorithmes  de  projection  rapides  afm 
d’estimer  le  sous-espace  des  brouilleurs  et  ainsi  calculer  I’estimateur  spectral.  On  decrit  trois 
algorithmes  le  premier  6tant  plus  efficace  du  point  de  vue  du  nombre  de  calculs  mais  moins 
performant  du  point  de  vue  du  pouvoir  de  resolution,  le  dernier  etant  moins  rapide  mais  plus 
performant.  Ces  algorithmes  donnent  un  meilleur  compromis  entre  le  rendement  et  la  charge  de 
calculs  que  les  algorithmes  connus  jusqu’i  present.  On  a  compare  le  rendement  de  ces  trois 
techniques  avec  celui  de  la  technique  MUSIC  pour  deux  scenarios  differents  qui  6valuent  le 
pouvoir  de  resolution  et  I’effet  de  la  puissance  des  brouilleurs. 
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EXECUTIVE  SUMMARY 


Accurate  angular  location  and  resolution  of  jamming  sources  is  important  for  military 
sensor  systems  such  as  ESM  and  radar.  This  report  describes  three  techniques  for  processing  the 
outputs  of  antenna  arrays  to  produce  improved  performance.  Specifically,  this  report  describes 
and  evaluates  computationally  efficient  algorithms  for  direction  finding;  these  reduce  the 
requirement  for  high  speed  computation  and  hence  the  cost  of  implementation. 

The  algorithms  to  be  described  and  evaluated  are  based  on  projection  techniques;  the 
jamming  signal  subspace  is  estimated  from  the  data  vectors  obtained  from  the  antenna  array 
elements.  The  jamming  signal  subspace  is  then  used  to  calculate  the  directional  spectrum,  which 
basically  consists  of  calculating  the  power  of  the  projection  of  the  search  vector  into  the  subspace 
orthogonal  to  the  jamming  signal  subspace.  As  the  search  vector  gets  closer  to  the  jammer 
direction,  the  power  of  the  projected  vector  decreases.  By  taking  the  inverse  of  the  projeC'Cd 
power  we  obtain  a  directional  spectrum  with  the  peak  positions  reflecting  the  jammer  directions. 

Three  algorithms,  all  based  on  data  selection  criteria  followed  by  Gram-Schmidt 
orthogonalization,  are  described  and  evaluated.  The  algorithms  are  named  as  follows:  (1)  DVO 
for  data  vector  orthogonalization,  (2)  DVSO  for  data  vector  selection  and  orthogonalization,  (3) 
DVSO-COVAR  for  repeated  applications  of  DVSO  followed  by  a  combining  of  the  selected  data 
vectors  into  a  covariance  matrix  which  in  turn  is  followed  by  a  final  application  of  DVSO  to  the 
columns  of  the  covariance  matrix. 

The  above  three  algorithms  are  evaluated  by  means  of  simulations  and  compared  with 
the  standard  MUSIC  technique  (MUSIC  stands  for  Multiple  Signal  Classification).  The  latter  is 
very  computationally  demanding.  In  terms  of  performance  versus  computational  load,  all  three 
techniques  compare  favourably  with  the  MUSIC  technique.  The  three  techniques  trade  off 
performance  against  computational  load,  DVO  being  the  fastest  in  terms  of  computations  but 
providing  poorest  resolution  capability  while  DVSO-COVAR  provides  best  resolution  capability 
but  at  a  cost  of  a  considerable  computational  load. 
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DATA  SELECTION  FOR  FAST  PROJECTION  TECHNIQUES: 

A  COMPARATIVE  STUDY  OF  DIRECTION  FINDING  PERFORMANCE 

I.O  INTRODUCTION 


Mcxlem  signal  processing  techniques  provide  the  capability  for  very  high-resolution 
determination  of  angles-of-arrival.  While  the  performance  of  these  techniques  is  very  good,  they 
have  very  high  computational  loads,  precluding  their  use  in  many  circumstances.  Projection 
techniques  based  on  eigenanalysis  form  the  basis  of  the  best  performing  and  most  popular  high- 
resolution  array  signal-processing  techniques.  Such  techniques  have  been  used  for  argle-of- 
arrival  estimation,  for  example  the  MUSIC  technique  of  Schmidt  [1]  and  the  equivalent 
formulations  of  Bienvenu  [2],  as  well  as  for  adaptive  antenna  nulling  of  jammers  [3,4].  Fast 
projection  techniques,  based  on  direct  orthogonalization  of  the  data  vectors,  avoid  time 
consuming  eigenanalysis  but  with -some  necessary  loss  of  performance  [5, 6, 7, 8].  Toulgoat  and 
Turner  [9,10]  have  introduced  aspects  of  data-vector  selection  to  obtain  a  better  trade-off 
between  performance  and  speed  for  adaptive  nulling  of  jamming  signals. 

This  report  is  an  evaluation  of  the  fast  projection  methods  of  Toulgoat  and  Turner  as 
applied  to  high-resolution  direction  finding  similar  to  the  evaluation  carried  out  for  adaptive 
nulling  in  [10],  As  in  [10],  we  consider  a  family  of  techniques  starting  with  the  fastest  but 
poorest  performing  and  progressing  to  the  slowest  but  best  performing.  We  consider  three 
techniques:  DVO  for  Data  Vector  Orthogonalization;  DVSO  for  Data  Vector  Selection  and 
Orthogonalization;  and  DVSO-COVAR,  where  DVSO  is  first  used,  the  selected  vectors  are  then 
formed  into  a  covariance  matrix,  and  DVSO  is  then  applied  to  the  columns  of  the  covariance 
matrix.  The  data  vectors  selection  technique  used  in  DVSO-COVAR  is  similar  to  a  technique 
of  Reilly  and  Law  [11].  However,  [11]  applies  the  selection  process  to  the  permutation  and 
selection  of  columns  of  a  covariance  matrix.  In  contrast,  we  apply  the  data  selection  prior  to  the 
calculation  of  the  covariance  matrix  in  order  to  choose  the  data  vectors  used  in  the  formation 
of  the  covariance  matrix  and  then  apply  a  procedure  which  is  equivalent  to  that  of  [11]  for 
selecting  a  subset  of  the  columns  vectors  of  the  covariance  matrix.  The  major  contributions  of 
this  report  are  the  new  algorithms  and  the  evaluation  of  their  performance  as  a  function  of 
computational  load  for  direction  finding. 

The  basis  of  the  evaluation  is  the  resolution  performance  achieved  versus  the 
computational  load  as  determined  using  Monte-Carlo  simulations.  In  order  to  have  a  consistent 
measure  of  resolution  performance,  we  define  two  signals  as  being  resolved  when  the  estimator 
has  two  peaks  indicating  the  two  signals  near  the  true  signal  positions  and  where  there  is  a  dip 
of  at  least  3  dB  between  the  peaks.  The  minimum  resolution  threshold  is  then  the  angular 
separation  between  the  signals  when  this  criterion  is  just  satisfied.  The  evaluation  is  carried  out 
for  two  test  cases:  T,  where  the  resolution  threshold  is  determined  for  fixed  and  high  signal-to- 
noise  ratio;  and  Tj  where  the  effect  of  different  signal-to-noise  ratios  on  the  directional  spectrum 
is  examined.  The  three  techniques  are  compared  among  themselves  as  to  their  performance 
versus  computational  loading.  The  final  comparison  is  then  made  with  the  MUSIC  algorithm. 
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2.0  SIGNAL  MODEL 


In  this  report,  we  consider  only  a  uniform  linear  array  with  equally  spared  elements. 
However,  the  techniques  developed  are  also  applicable  to  arbitrarily  shaped  arrays  of  non- 
uniformly  spaced  elements.  The  signals  of  interest  comprise  jamming  interference  received  by 
the  radar  system.  Such  signal  sources  can  be  considered  to  be  completely  uncorrelated.  The 
objective  is  to  resolve  these  sources  and  to  accurately  determine  their  directions. 

Radar  systems  are  usually  narrow  band;  this  means  that  the  received  signal  wave  front 
has  nearly  perfect  spatial  correlation  over  the  face  of  the  array,  provided  that  the  time  taken  to 
cross  the  array  is  small  compared  with  the  reciprocal  of  the  system  bandwidth.  The  n*^  data 
vector  at  time  t„  is  expressed  as  : 

L 

i-1 


where  L  is  the  number  of  signals,  n„  is  a  receiver-noise  vector  with  mutually  independent 
components  and  power,  E{  |  n„  p}  =  K  with  K  the  number  of  array  elements.  The  quantity 
ji(tj  is  a  complex  Gaussian  random  variable  representing  the  i"*  signal  amplitude  at  time  t„,  and 
Si  is  a  deterministic  vector  representing  the  direction  of  arrival  of  the  i"*  signal  defined  as 

=  [1,  exp  (j2itdsin(0^)  /X) ,  .  .  . ,  exp  {j2iid{K-l)sind^/X)  ]  ^ 

Here  d  is  the  inter-element  spacing,  and  0;  is  the  direction  of  arrival  of  the  i*  signal.  The  sample 
vectors  or  "snapshots"  are  taken  at  time  intervals  such  that  jjftj  and  ji(t„.i)  are  completely 
independent. 


3.0  PERFORMANCE  EVALUATION 
3.1  MONTE-CARLO  SIMULATIONS 

Monte-Carlo  simulations  are  used  to  evaluate  the  performance  of  the  various  techniques. 
One  hundred  independent  trials  are  carried  out  for  each  parameter  setting  under  test  with 
performance  measures  averaged  over  these  trials.  The  jamming  signals  are  generated  as  samples 
of  a  zero-mean  complex  Gaussian  variate  having  perfect  spatial  correlation  over  the  array  and 
zero  temporal  correlation  from  one  data  vector  to  the  next.  This  gives  rise  to  a  signal  amplitude 
that  is  Rayleigh  distributed.  Receiver  noise  is  also  generated  as  samples  of  a  zero-mean  complex 
Gaussian  variate  but  having  zero  spatial  correlation  over  the  array  and  zero  temporal  correlation 
from  one  data  vector  to  the  next.  The  receiver  noise  power  at  the  element  level  is  set  at  = 
10^. 


The  results  are  presented  for  ten-element  and  forty-element  arrays  of  omnidirectional 
elements  with  an  element  spacing  of  one  half  wavelength. 
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3.2  PERFORMANCE  MEASURES 


3.2.1  DIRECTIONAL  SPECTRUM 


The  directional  spectrum  F(6)  used  for  angle-of-arrival  estimation  is  given  as 


F(0)  =  - - - 

fl(0)"P  3(6) 


(1) 


where  P  is  the  projection  matrix 

(2) 


s(0)  is  defined  as 

a  (6)  =  [  1 .  exp  ( s in0 )  .... ,  exp(  — ^sin0)  ]  ^ 

d  is  the  spacing  between  array  elements,  X  is  the  signal  wavelength  and  d  is  the  angle  searched. 

The  entire  space  is  searched,  using  low-resolution  Fourier  beamforming  prior  to  the 
application  of  high-resolution  spectral  estimation  to  localize  the  regions  of  interest.  This 
procedure  is  used  to  reduce  the  number  of  directions  (denoted  by  nj  that  are  investigated.  The 
present  state  of  the  art  permits  a  resolution  of  about  V*  of  a  beamwidth.  The  minimum  number 
of  directions  to  be  searched  is  therefore  n4=4  directions  within  a  beamwidth.  In  practice 
however,  the  number  may  not  be  large  enough  to  resolve  the  peak  positions.  In  this  report,  we 
consider  nd=23.  This  value  is  larger  than  what  might  be  needed  in  practice;  however,  such  a 
large  value  of  suits  our  purpose:  to  show  that  the  fast  projection  algorithms  can  be  used  for 
direction  flnding  and  to  compare  their  resolution  with  that  of  the  MUSIC  method. 


3.2.2  MINIMUM  RESOLUTION  THRESHOLD 

The  performance  measure  used  is  the  minimum  resolution  threshold,  A6^,  for  two 
uncorrelated  stochastic  signals  having  the  same  Signal-to-Noise  Ratio  (SNR).  Here  the  term  SNR 
indicates  the  ratio  of  the  power  of  the  signal  to  the  power  of  the  receiver  noise  measured  at  the 
element  receiver  output.  In  each  simulation  trial,  the  two  sources  were  considered  resolved  if 
they  met  the  two  following  requirements: 
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a)  The  spectral  output  of  the  algorithm  under  test  had  two  peaks  at  or  near  the  true  source 
locations.  An  error  of  ±(02-0i)/lO  in  location  was  allowed. 

b)  The  peaks  at  and  62  were  more  than  3  dB  above  the  value  of  ¥(d)  at  6=6,  ±  A6/2, 
where  i=l,2  and  A6=62-6i. 

This  measure  of  angular  resolution  is  given  in  beamwidths  (BW).  One  beam  width  is 
defined  as  X/D  where  X  is  the  signal  wavelength  and  D  is  the  antenna  aperture. 


3.2.3  ERROR  IN  THE  ESTIMATION  OF  DIRECTIONS 

Since  the  data  vectors  contain  receiver  noise  in  addition  to  signals,  the  estimate  of  the 
signal  subspace  is  not  perfect.  This  causes  error  in  the  determination  of  signals  directions.  As 
seen  in  3.2.2,  an  error  of  ±A0/1O  was  allowed  in  the  calculation  of  the  minimum  resolution 
threshold.  To  determine  this  error  more  precisely,  we  calculate  the  value  of  F(0)  at  intervals  of 
A0/1(X)  to  obtain  a  precision  of  1  %  of  A6.  If,  for  example,  we  observe  a  peak  at  0,  +  3A0/1OO, 
we  conclude  that  the  position  is  off  by  2.5  to  3.5%  of  the  angle  separation,  A6. 


3.3  TEST  CASES 

We  have  specified  two  test  cases,  T,  and  Tj,  which  are  designed  to  determine  the 
minimum  threshold  resolution,  A0„i„,  achieved  by  the  various  techniques  and  to  test  the  effect 
of  the  signal  strength  on  the  directional  spectrum,  respectively. 


3.3.1  TEST  CASE  T, 

Here  we  have  specified  the  Signal-to-Noise  Ratio,  SNR,  at  the  receiving  element  output 
to  be  40  dB.  The  number  of  signals  is  designated  as  L=2.  The  first  signal  is  set  at  0,  =20°  while 
the  second  signal  direction  is  at  62=6^+A6.  The  separation  A6  is  then  reduced  until  the  threshold 
resolution,  A6^,  is  reached. 


3.3.2  TEST  CASE  Tj 


Here  we  have  defined  two  equi-powered  signals  at  0,  =20°  and  02=30°.  We  consider 
four  values  of  SNR  of  10,  20,  30  and  40  dB  as  measured  at  the  element  level.  This  test  case  is 
used  to  examine  the  performance  of  the  various  estimators  as  a  function  of  SNR. 
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4.0  DATA  VECTOR  ORTHOGONALIZATION 


The  Data  Vector  Orthogonalization  (DVO)  technique  is  an  efficient  method  for 
constructing  a  basis  for  the  (jamming)  signal  subspace  and  computing  the  corresponding 
directional  spectrum  F(6).  A  set  of  M  basis  vectors,  (v„},  is  generated  directly  from  a  set  of 
sample  vectors,  {y„:  m  =  l  to  M},  selected  from  a  larger  group  of  N  data  vectors,  {x„:  n  =  l  to 
N},  using  the  Gram-Schmidt  (GS)  orthonormalization  method.  Two  thresholds  are  included  in 
the  procedure:  an  internal  threshold  A  that  rejects  data  vectors  which  provide  little  additional 
information  about  the  jamming  signal  subspace  and  an  external  threshold,  Aj,  which  forces  the 
continued  collection  of  data  vectors  and  provides  a  termination  criterion  for  the  algorithm.  The 
use  of  an  external  threshold  in  this  manner  was  proposed  by  Nickel  [7]. 

The  m*  basis  vector  v„  is  generated  from  the  n*  data  vector  x„  using 

/n-1 

-  E 

J-1 


We  then  compare  the  squared  magnitude  of  ii„  to  an  internal  threshold  given  by  A=2K<r^.  The 
objective  is  to  eliminate  data  vectors  whose  information  content  is  too  low.  If  ju„,i^  >  A,  then 
u„  is  accepted  and  a  new  basis  vector,  v„,  is  generated  as  v„  =  u„/ 1  u„  | .  If  is  not  accepted, 
the  procedure  is  started  again  with  a  new  data  vector  x„. 

Following  the  generation  of  v„,  a  statistic  T  is  calculated  as  in  [7] 

1 

n|y/ 

i-1 

and  tested  against  an  external  threshold,  Aj.  Matrix  Y  is  constructed  with  the  m  data  vectors 
selected  from  {x„}: 

When  T  '  Aj,  the  signal  subspace  has  been  sufficiently  well  estimated  and  the  procedure  is 
stopped.  Otherwise,  the  procedure  continues  with  a  new  data  vector,  x„.  The  accuracy  of  the 
estimation  of  the  signal  subspace  depends  on  the  choice  of  A^.  Indeed,  there  is  a  trade-off 
between  the  accuracy  of  the  estimate  and  the  computational  load  depending  on  the  choice  of  A^. 
The  number  of  data  vectors,  N,  that  are  examined  and  the  number,  M,  selected  are  random 
variables  for  a  given  threshold  value.  At.  The  dependence  of  the  computational  load  on  N  and 
M  will  be  evaluated  by  means  of  Monte-Carlo  simulations. 
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4.1  COMPUTATIONAL  LOAD 


The  computational  count  comprises  two  terms;  the  first  term,  C^,  is  associated  with  the 
calculation  of  the  basis  vectors  {v„};  the  second  term,  C,*,  represents  the  number  of  complex 
multiplications  required  for  the  calculation  of  the  directional  spectrum  F(0). 

In  the  DVO  technique,  we  assume  that  N  data  vectors  are  used  to  compute  M 
orthonormal  vectors  spanning  the  signal  subspace.  An  upper  bound  to  the  number  of  complex 
multiplications  required  is  then 

^DVO  ~  ■*"  ^se 


where 


Cfc  =  KN^ 


0  .  SKN  +  0  .  SKM  + 


(N-1)  iK*l) 
2 


and 


=  nj  M(K+0.5)  (3) 

where  is  the  number  of  angles  searched  i.e.  the  number  of  points  at  which  the  directional 
spectrum,  F(5),  is  evaluated. 


4.2  COMPUTER  SIMULATIONS 

We  first  determine  the  minimum  resolution  threshold,  ^6^,  of  the  DVO  algorithm  for 
arrays  of  10  and  40  elements  using  test  case  T,.  Figure  1  shows  the  resolution  threshold,  A6^, 
versus  the  external  threshold  Aj.  As  Aj  increases,  the  accuracy  of  the  estimate  of  the  signal 
subspace  decreases  and  therefore  Ad„i„  increases.  The  errors  in  the  estimation  of  the  directions 
of  the  sources  vary  between  0%  and  6.5%  of  A6  (Aff=A0^  in  this  case)  for  both  sources  and 
both  arrays. 

Next,  we  analyze  the  effect  of  the  signal  strength.  Figure  2  shows  the  directional 
spectrum  F(0)  versus  $  for  a  threshold  =  0-5  using  test  case  Tj.  The  peak  heights  are 
proportional  to  the  relative  SNR  for  each  signal.  The  error  in  estimating  the  directions  is  less 
than  0.5%  of  AO  for  SNR  of  20,  30  and  40  dB  and  is  less  than  1.5%  of  AO  for  SNR =10  dB. 
Further  computer  simulations  showed  that  as  the  signal  strength  increases,  the  minimum 
resolution  threshold,  AO^,  achieved,  becomes  lower.  As  an  example,  two  sources  at  =  20° 
and  02  =  25  °  are  resolved  with  a  threshold  Ax =0.5  for  SNR =40  dB  but  cannot  be  resolved 
even  with  a  threshold  Ax=lxl0'*  if  SNR=10  dB.  In  the  latter  case,  it  was  not  possible  to  find 
an  estimate  of  the  signal  subspace  sufficiently  accurate.  As  expected,  the  accuracy  of  the 
estimate  of  the  signal  directions  increases  with  SNR. 
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Fig.  1  versus  for  DVO  with  L=2  and 

SNR=40  dB  for  arrays  of  10  and  40  elements. 
Monte-Carlo  simulations  with  100  trials. 
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Fig.  2  F(0)  versus  0  for  DVO  with  K=10, 

L=2,  Aj=0.5.  Monte-Carlo  simulations  with  100 
trials. 


5.0  DVSO  TECHNIQUE 

The  DVSO  technique  is  based  on  a  selection  of  the  best  M  vectors  from  a  larger  group 
of  N  vectors  which  are  sampled  when  only  the  jamming  signals  are  present  and  stored  before 
the  process  begins.  The  selected  M  vectors  are  best  in  the  sense  that  they  yield  the  most  accurate 
estimate  of  the  signal  subspace  as  compared  with  any  other  set  of  M  vectors  selected  from  the 
larger  set  of  N.  In  the  algorithm  for  estimating  the  signal  subspace,  we  use  the  following 
notation:  {x„}  comprises  the  set  of  N  sample  vectors,  {u„}  denotes  the  M  best  sample  vectors 
chosen  from  the  {x„},  and  {v„}  denotes  the  orthonormaJ  basis  formed  from  the  {u,„}  by  the  GS 
method. 

The  data  selection  is  integrated  with  the  GS  orthonormalization  method.  The  data  vector 
chosen  has  the  largest  projected  amplitude  in  the  subspace  orthogonal  to  the  previously  selected 
(m-1)  basis  vectors  {v,,  ...,  v„}.  The  squared  magnitude  of  this  projection  is 

<  =  n-1 . N 

where  is  the  projection  matrix  determined  from  the  basis  {v,,  ...,  as  in  equation  (2). 
The  iterations  are  initiated  with  Po=I,  the  identity  matrix. 
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The  m*  basis  vector  v„  is  generated  from  the  n*  data  vector  {x„}  using 

“m  =  -  E  'O- 

J-i 


We  then  compare  the  squared  magnitude  of  u„  to  a  threshold  given  by  A=2Ka^.  If  ju„j^  >  A 
then  v„=ii„/  j  Un,|  and  the  process  starts  again  with  another  data  vector.  Otherwise,  the  algorithm 
terminates  with  an  estimate  of  the  signal  subspace  given  by  where  M  equals  the  final 

value  of  m  in  the  iteration  minus  1.  The  angle  estimator  is  then  obtained  using  equation  (1). 


5.1  COMPUTATIONAL  COUNT 

Although  the  concept  of  the  projection  matrix  is  central  to  the  derivation  of  the  technique, 
it  is  never  explicitly  calculated.  It  is  computationally  more  efficient  to  work  directly  with  the 
basis  vectors.  The  projection  matrix  is  only  required  in  the  selection  of  the  data  vector.  We  use 
the  following  properties  of  projection  matrices  to  simplify  the  expression  for 

P  “  =  p 
and 

p  .p  .  =  p  . 


We  can  now  write 


e 


m 

n 


m-l 


Jc-l 


We  obtain  the  iteration  for  as 


With  those  simplifications,  the  computational  load  for  the  DVSO  technique  is  given  by: 
CdvSO  “  C|,  +  C„ 

where 

Cfc  =  K  (0.5  +  0.5N  +  2M  +MN  +  M^)  +  0.5  MN 
and  C^  is  given  by  equation  (3) 
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5.2  COMPUTER  SIMULATIONS 


We  first  determine  the  minimum  resolution  threshold,  A6^,  for  the  DVSO  technique 
using  test  case  T,.  Figure  3  shows  Ad^  versus  N  for  arrays  of  10  and  40  elements.  As  N 
increases,  the  resolution  threshold  decreases  down  to  a  limit  which  is  around  0.12  BW  for 
K=10  and  0.13  BW  for  K=40.  Increasing  N  would  slightly  lower  the  minimum  resolution 
threshold.  N  has  been  limited  to  20  in  the  results  presented  here.  It  is  observed  that  as  K 
increases,  the  resolution  in  BW  increases  by  a  small  amount.  We  then  consider  the  error  on  the 
positions  of  the  directions  for  the  results  of  Figure  3.  The  error  in  estimating  the  source 
directions  lies  between  15%  tn  4  5%  of  A6  for  0,  and  2.5%  to  5.5%  of  A6  for  62  for  both 
arrays,  for  Ad  near  the  resolution  limit,  Ad^. 

We  next  analyze  the  effect  of  signal  strength  on  the  directional  spectrum.  Figure  4  shows 
the  directional  spectrum  F(fl)  versus  d  using  test  case  Tj  with  N=5  for  SNR=  10,  20,  30,  40  dB. 
The  peaks  of  the  function  reflect  the  signal  strength.  Other  computer  simulations  (results  not 
given)  showed  that  the  height  of  the  peaks  is  fairly  insensitive  to  N,  provided  N  is  high  enough 
to  resolve  the  sources.  It  is  also  observed  that  the  resolution  threshold  depends  on  the  signal 
strength.  As  an  example,  for  ^,=20“  and  ^2=25°  we  cannot  find  an  N<20  for  which  the 
sources  are  resolved  for  SNR  =  10  dB  but  they  are  resolved  with  N=3  for  SNR=40  dB. 

The  error  in  estimating  the  signal  directions  has  been  determined  for  test  case  Tj. 
The  results  showed  that  the  error  is  less  than  1.5%  of  Ad  for  all  SNR’s.  We  observe  from  test 
cases  T,  and  T2  that,  as  the  angle  separation,  Ad,  increases,  the  error  in  the  signal  direction 
decreases. 


6.0  COMPARISON  OF  DVO  AND  DVSO 

Both  the  DVO  and  DVSO  methods  combine  GS  orthogonalization  with  a  data  selection 
process  for  the  computation  of  the  directional  spectrum  F(0).  In  the  DVO  method,  the  data 
selection  process  consists  of  eliminating  those  data  vectors  which  provide  little  additional 
information  about  the  jamming  signal  subspac&  while  continuing  to  take  new  sample  vectors  until 
a  second  external  threshold  criterion  is  satisfied.  The  data  selection  is  carried  out  after  the  GS 
orthogonalization.  In  the  DVSO  technique,  the  data  selection  process  consists  of  choosing  the 
data  vector  which  has  the  largest  projected  amplitude  into  the  subspace  orthogonal  to  the 
previously  selected  (m-l)  basis  vectors  {v„  ...,  The  data  selection  process  of  the  DVO 
method  is  not  as  efficient  as  that  of  the  DVSO  method.  In  the  DVO  method,  we  eliminate  data 
vectors  that  are  not  good  enough  while  in  DVSO  we  select  the  best  data  vectors.  Therefore,  the 
DVSO  method  gives  a  more  accurate  estimate  of  the  signal  subspace  than  does  the  DVO 
technique;  the  data  selection  criterion,  however,  requires  more  computations. 

The  curves  of  performance  versus  number  of  complex  multiplications  are  obtained  as 
follows:  (1)  for  DVO  we  decrease  the  threshold,  Aj,  which  increases  the  number  of  calculations 
but  gives  better  performance;  (2)  for  DVSO  we  increase  N,  increasing  the  computational  load 
and  improving  the  resolution. 
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Fig.  3  versus  N  for  DVSO  with  L=2, 

SNR=40  dB.  Monte-Carlo  simulations  with  100 
trials. 


Fig.  4  F(0)  versus  0  for  DVSO  with  K=10, 

L=2  and  N=5.  Monte-Carlo  simulations  with  100 
trials. 


Figure  5  shows  the  minimum  resolution  threshold,  versus  complex  multiplications 
for  both  DVO  and  DVSO  techniques  for  an  array  of  10  elements.  Two  regions  of  operation  are 
observed.  Below  152  BW,  the  DVSO  method  performs  best  in  the  sense  that  it  requires 

less  computations  to  provide  a  given  Above  152  BW,  the  DVO  method  performs 

best.  These  results  are  obtained  with  nj=23  directions.  Figure  6  gives  the  same  performance  for 
an  array  of  40  elements.  We  notice  that  there  is  no  difference  between  DVO  and  DVSO  for  this 
array.  Since  the  computational  load  is  composed  of  two  terms  where  only  one  term  depends  on 
the  parameter  n^,  it  is  interesting  to  plot  Aff  against  the  two  terms  separately.  Figures  7  and  8 
shows  A0^  versus  Cb  and  0,^,  respectively,  for  an  array  of  40  elements.  From  Figure  7,  we  see 
that  the  calculation  of  the  basis  by  DVSO  is  more  computationally  intensive  than  that  of  the 
DVO  method.  Figure  8  shows  that  the  search  operation  is  less  computationally  intensive  for 
DVSO  when  Ad„^  is  less  than  0. 17  BW  while  requiring  about  the  same  number  of  computations 
as  DVO  for  A0^  greater  than  0. 17  BW.  Considering  both  terms  of  the  computational  load,  the 
DVSO  technique  is  advantageous  for  small  values  of  A0^  and  small  arrays  while  both  performs 
approximately  the  same  as  DVO  for  larger  Aff  and  larger  arrays. 
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Fig.  6  A0„,„  versus  complex  multiplications  for 

DVO  and  DVSO  with  K=40.  L=2  and  SNR=40 
dB.  Monte-Carlo  simulations  with  100  trials. 
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Fig.  7  A0„^,  versus  Q  for  DVO  and  DVSO 

with  K=40,  L=2  and  SNR=40  dB.  Monte-Carlo 
simulations  with  100  trials. 
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Fig.  8  A0„i„  versus  C„  for  DVO  and  DVSO 

with  K=40,  L=2  and  SNR=40  dB.  Monte-Carlo 
simulations  with  100  trials. 


7.0  DVSO  FOR  ESTIMATING  THE  COVARIANCE  MATRIX 


In  this  technique  we  employ  the  DVSO  method  twice.  The  first  step  is  to  subdivide  a  set 
of  QN  data  vectors  into  Q  groups  of  N  vectors.  DVSO  is  then  applied  to  select  the  best  M,  data 
vectors  out  of  the  i*  group,  i  =  1  to  Q.  This  is  repeated  Q  times  to  give 

0 

"r=  E"i 

i-i 


data  vectors  which  are  used  to  calculate  the  covariance  matrix  R; 


Wj. 


ri-i 


We  are  thus  able  to  obtain  an  estimate  of  the  covariance  matrix  using  fewer  but  higher 
quality  data  vectors  than  in  the  conventional  approach. 

The  final  step  applies  the  DVSO  method  to  the  columns  of  the  covariance  matrix  to  find 
an  orthonormal  basis  and  calculate  a  directional  spectrum.  This  method  is  called  DVSO- 
COVAR. 


The  data  vector  selection  technique  used  is  similar  to  a  technique  of  Reilly  and  Law  [11]. 
However,  [11]  applies  the  selection  process  to  the  permutation  and  selection  of  columns  of  a 
covariance  matrix.  In  contrast  we  apply  the  criterion  to  the  data  vectors  themselves  thus  reducing 
the  number  of  data  vectors  used  for  computing  the  covariance  matrix  and  therefore  reducing  the 
computational  load.  When  we  calculate  the  covariance  matrix  as  in  the  DVSO-COVAR 
algorithm,  the  data  vectors  are  subjected  to  the  selection  process  before  being  combined  in  the 
covariance  matrix.  We  then  apply  a  procedure  which  is  equivalent  to  that  of  [1 1]  for  selecting 
a  subset  of  the  columns  vectors  which  are  used  to  generate  a  basis  for  the  signal  space  by  means 
of  the  Gram-Schmidt  process. 

The  number  of  computations  required  by  the  DVSO-COVAR  technique  is  data-dependent; 
the  mean  number  of  calculations  is  estimated  by  Monte-Carlo  simulations.  For  a  single  trial 
evaluation  of  the  directional  spectrum,  the  number  of  multiplications  is 

i-l  ^ 


where  C^  is  the  number  of  computations  for  the  calculation  of  the  basis  by  the  DVSO  technique 
and  C^  is  the  number  of  computations  for  the  evaluation  of  the  directional  spectrum;  H  is  the 
number  of  basis  vectors  used  to  calculate  the  directional  spectrum. 
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7.1  COMPUTER  SIMULATIONS 


We  first  analyze  the  effect  of  varying  the  parameters  Q  and  N  on  the  minimum  resolution 
threshold  A6^.  Figure  9  shows  the  minimum  resolution  threshold,  Ad^,  versus  N  for  an  array 
of  10  elements  with  Q=  1,2,3  and  4  using  test  case  T,.  As  for  DVSO,  increasing  N  lowers  the 
minimum  resolution  threshold.  This  threshold  is  also  lowered  by  increasing  Q.  For  Q=l,  the 
minimum  resolution  threshold  is  the  same  as  that  obtained  for  the  DVSO  technique  while  for 
Q=4,  Ad^  can  be  as  low  as  0.09  BW.  As  in  the  DVSO  technique,  N  is  limited  to  20.  Figure 
10  shows  the  minimum  resolution  threshold,  Ad^,  versus  complex  multiplications.  We  see  that 
increasing  Q  decreases  the  minimum  resolution  threshold  achieved  but  for  a  given  Ad^,  the 
computational  load  is  about  the  same  for  the  different  values  of  Q.  This  can  be  explained  by  the 
fact  that  as  Q  increases,  we  need  a  lower  N  to  achieve  the  same  resolution.  Therefore  for  small 
values  of  A6^,  there  is  no  difference  in  the  computational  load.  The  error  in  estimating  the 
directions  of  the  sources  for  the  results  of  Figure  10  (test  case  T,)  lies  between  0.5%  to  6.5% 
of  Ad  for  6i  and  1.5%  to  5.5%  of  Ad  for  (here  Ad=Ad^. 

Figure  1 1  shows  the  minimum  resolution  threshold,  Ad^  versus  Q  for  different  values 
of  N  using  test  case  T,  for  an  array  of  10  elements.  Increasing  N  for  a  given  Q  lowers  the 
minimum  resolution  threshold.  The  decrease  in  Ad^  obtained  by  increasing  N  is  very  small  (less 
than  0.01  BW).  Figure  12  shows  Ad^  versus  complex  multiplications  for  given  values  of  N  (the 
increase  in  the  number  of  multiplications  for  a  given  N  is  due  to  the  increase  in  the  value  of  Q). 
Up  to  a  certain  point,  increasing  N  does  not  dramatically  increase  the  computational  load  while 
permitting  a  better  resolution.  However,  the  curve  of  resolution  versus  calculation  load  quickly 
flattens  out,  so  that  a  very  high  value  of  N  does  not  provide  a  concomitant  improvement  in  the 
resolution  threshold.  The  bias  of  the  estimates  of  signal  directions  varies  from  1.5%  to  6.5% 
of  Ad  for  estimates  of  A,  and  from  2.5%  to  6.5%  of  Ad  for  estimates  of  02- 

We  now  compare  the  results  obtained  for  K=  10  with  those  obtained  for  an  array  of  40 
elements.  Figure  13  shows  Ad^  in  BW  versus  N  for  DVSO-COVAR,  Q=4  with  K=10  and 
K=40  using  test  case  T,.  The  curves  are  quite  similar  even  though  the  resolution  in  BW  is 
slightly  better  for  a  smaller  array. 

Similarly,  we  plotted  Ad^  versus  Q  for  DVSOCOVAR,  N=5  for  arrays  of  10  and  40 
elements  in  Figure  14.  The  two  curves  lie  on  the  top  of  each  other. 

We  next  analyze  the  effect  of  signal  strength  on  the  directional  spectrum,  F(fl),  using  test 
case  T2.  Figure  15  shows  F(fl)  versus  d  for  Q=  1  and  different  values  of  SNR.  The  peaks  reflect 
the  signal  strength  with  an  offset  of  5  dB.  Figure  16  shows  the  directional  spectrum  F(d)  versus 
d  for  SNR =40  dB  and  for  different  values  of  Q.  As  Q  increases,  the  peaks  become  higher. 

The  position  error  is  less  than  0.5%  of  Ad  for  every  value  of  SNR  and  Q.  Moreover,  the 
sharpness  of  the  peaks  increases  with  Q  for  a  given  SNR.  Finally,  the  sharpness  of  the  peaks 
increases  with  SNR  for  a  given  Q.  As  for  DVSO,  we  conclude  from  test  cases  Tj  and  T2  that 
the  error  in  the  estimate  of  the  signal  directions  decreases  with  increasing  Ad. 
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Fig.  9  versus  N  for  DVSO-COVAR  with 
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Fig.  1 1  versus  Q  for  DVSO-COVAR  for 
N=S,  10  and  20  with  K=10,  L=2  and  SNR=40 
dB.  Monte-Carlo  simulations  with  100  trials. 
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Fig.  12  versus  complex  multiplications  for 
DVSO-COVAR  with  K=I0.  L=2  and  SNR=40 
dB.  Monte-Carlo  simulauins  with  100  trials. 


8.0  COMPARISON  OF  DVSO  AND  DVSO-COVAR 


We  first  compare  the  DVSO  and  DVSO-COVAR  methods  for  test  case  Tj  on  the  basis 
of  their  resolution  versus  computational  load.  Figure  17  shows  Ad^  versus  complex 
multiplications  for  DVSO  and  DVSO-COVAR  with  nj=23  search  directions  for  K=10.  Two 
curves  are  shown  for  DVSO-COVAR:  one  for  a  fixed  Q  (Q=4)  with  N  varying,  the  other  for 
fixed  N  (N=5)  with  Q  varying.  There  are  two  regions  of  operation:  for  Ad^  >  0.124  BW, 
DVSO  performs  best;  for  Ad^  <  0.124  BW,  DVSO-COVAR  performs  best.  From  Figure  17 
we  see  that  DVSO-COVAR  achieves  a  much  lower  resolution  threshold  by  varying  Q  for  a  fixed 
N  than  by  varying  N  for  a  fixed  value  of  Q.  Figure  18  shows  Ad^  versus  the  number  of 
complex  multiplications  for  both  techniques  with  nj=250.  As  n^  increases,  DVSO-COVAR 
becomes  slightly  more  efficient  for  low  values  of  A8^.  For  high  values  of  Ad^  where  the  value 
of  N  is  at  its  smallest,  the  advantage  of  DVSO  is  evident.  It  should  be  noted  however,  that 
DVSO-COVAR  can  achieve  a  very  high  resolution  -  a  resolution  not  achievable  with  DVSO- 
albeit  at  the  cost  of  a  considerable  computational  load.  We  next  analyze  the  contribution  of 
the  terms  C^  and  to  the  computational  load.  Figure  19  shows  A6^  versus  the  number  of 
complex  multiplications  required  for  the  evaluation  of  the  basis,  C^,  for  both  techniques  in  the 
K  =  10  case.  As  expected,  the  DVSO  method  is  more  computationally  effective  than  DVSO- 
COVAR  for  the  determination  of  the  basis.  Figure  20  shows  A6^  versus  the  number  of 
computations  required  for  the  evaluation  of  the  directional  spectrum,  C,,.  for  the  same  techniques 
with  n„=23.  Here  we  see  that  for  low  values  of  A6^,  DVSO-COVAR  requires  less 
computations  than  DVSO  but  for  high  values  (above  0.2  BW)  the  number  of  computations 
required  is  the  same. 

Finally,  Figure  21  shows  the  DVSO  and  DVSO-COVAR  methods  for  K=40  with  0^=23 
using  test  case  T,.  For  AO^  below  0,13  BW,  DVSO-COVAR  performs  better  while  for  Ad^ 
above  0.13  BW  the  DVSO  method  performs  better.  We  conclude  that  for  direction  finding,  it 
is  advantageous  to  use  the  DVSO  method  when  only  a  moderate  degree  of  superresolution  is 
required.  These  conclusions  are  valid  for  a  small  number  of  signals  such  as  L=2.  It  is  likely  that 
if  the  number  of  jammers,  L,  was  increased  in  comparison  to  K,  the  DVSO-COVAR  method 
would  give  better  results. 
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Fig.  17  versus  complex  multiplications  for 
DVSO  and  DVSO-COVAR  with  K=10.  L=2. 
n^23  and  SNR=40  dB.  Monte-Carlo  simulations 
with  100  trials. 


Fig.  18  A0„,„  versus  complex  multiplications  for 
DVSO  and  DVSO-COVAR  with  K=10.  L=2. 
SNR=40  dB  and  0,1=250.  Monte-Carlo 
simulations  with  100  trials. 
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Fig.  19  versus  C»  for  DVSO  and  DVSO- 
COVAR  with  K=10.  L=2,  SNR=40  dB.  Monte- 
Carlo  simulations  with  100  trials. 


Fig.  20  A0^  versus  for  DVSO  and  DVSO- 
COVAR  with  K=10,  L=2.  SNR=40  dB.  Montc- 
Carlo  simulations  with  100  trials. 


Complex  multiplications 


Fig.  21  versus  complex  multiplications  for 
DVSO  and  DVSO-CCVAR  Q=4  with  L=2  and 
SNR=40  dB.  Monte-carlo  simulations  with  100 
trials. 


9.0  COMPARISON  WITH  MUSIC 

Finally,  we  compare  the  fast  projection  techniques  with  the  MUSIC  method.  Since  we 
used  a  uniformly  spaced  linear  array,  the  root-MUSIC  technique  [12]  would  have  been  a  more 
appropriate  choice  but  for  the  sake  of  generality  we  used  the  standard  MUSIC  technique  because 
our  techniques  are  also  applicable  to  non-linear  arrays. 


9.1  DESCRBPnON  OF  MUSIC  TECHNIQUE 

In  the  MUSIC  algorithm  [1],  the  jamming  signal  subspace  is  constructed  with  the 
strongest  eigenvectors  of  the  covariance  matrix.  The  covariance  matrix  is  calculated  from  S  data 
vectors  as  follows: 


Jt  = 


where  Xj  is  the  i"*  data  vector  and  S  is  the  number  of  data  vectors  used  in  the  calculation  of  the 
covariance  matrix.  This  Hermitian  matrix  is  decomposed  in  terms  of  its  eigenvectors  as 
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where  X,  >  >  ...  >Xk  is  the  set  of  ordered  eigenvalues  and  {v,,...,Vk}  are  the  corresponding 

eigenvectors. 


The  Akaike  Information  Criterion  (AIC)  [13]  is  used  to  select  the  signal  space 
eigenvectors.  The  number  of  these  vectors  is  the  p  which  minimizes  the  function: 


AIC[p]  =  (JC-plln 
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9.2  COMPUTATIONAL  LOAD 

The  expression  for  the  computational  load  of  the  MUSIC  technique  comprises  three 
terms.  The  first  term,  Cc„,  is  the  number  of  computations  required  for  the  determination  of  the 
covariance  matrix: 

r  -  ■S'(XD  (X*l} 

'•can  ^ 


where  S  is  the  number  of  data  vectors  used  in  the  formation  of  the  covariance  matrix  and  K  is 
the  number  of  array  elements.  The  second  term,  C„,  is  associated  with  the  calculation  of  the 
eigenvalues  and  associated  eigenvectors.  A  very  efficient  way  to  compute  the  eigenvalues  of  a 
complex  Hermitian  matrix  (as  the  covariance  matrix)  employs  Householder  transformations  to 
reduce  the  matrix  to  an  equivalent  symmetric  tridiagonal  matrix.  The  implicit  QL  algorithm  [14] 
is  then  used  to  compute  the  eigenvalues  of  this  tridiagonal  matrix.  The  number  of  operations, 
Cu,  required  for  the  tridiagonalization  is  given  by 
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Since  the  average  number  of  iterations  per  eigenvalue  is  variable  in  the  QL  algorithm,  we  can 
only  provide  an  approximate  expression  for  the  computations.  The  average  number  of  iterations 
per  eigenvalue  is  typically  1.3-1. 6  [15].  We  chose  1.3  for  our  calculation  in  order  to  obtain  a 
lower  bound  on  the  computational  load  for  MUSIC.  Then  the  number  of  computations  for  the 
QL  algorithm  is 

-  6K^-3K 
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We  then  use  the  power  method  [15]  to  calculate  the  eigenvectors  associated  with  the 
strongest  eigenvalues.  The  number  of  operations  required  for  the  power  method,  Cp^,  is 

^pow  “  +  8nl^ 

where  n„  is  the  number  of  eigenvectors  to  be  calculated  and  f,(n„^)  is  the  number  of 
computations  required  for  the  full  eigen  solution  of  a  n^v  x  n^v  matrix.  We  assumed  that  the 
eigenvectors  are  obtained  after  six  iterations. 

The  total  number  of  computations  for  the  evaluation  of  the  eigenvalues  and  eigenvectors 
is 

^CT  ^pow 

The  last  term,  C.^,  is  associated  with  the  search  of  the  peaks  of  the  directional  spectrum: 

=  n^piK  +0.5) 


The  total  computational  load  is  then  given  by 

^MUSIC  ~  ^cm  ^0V  *  ^80 


In  the  next  section  we  will  show  that  the  projection  techniques  compare  favourably  to  the 
MUSIC  technique. 


9.3  SIMULATION  RESULTS 

We  compared  the  DVSO,  DVSO-COVAR  and  MUSIC  techniques  on  the  basis  of  test 
case  T|.  The  performance  of  the  MUSIC  technique  is  shown  in  Figure  22;  is  plotted  versus 
S  for  arrays  of  10  and  40  elements.  As  expected,  the  resolution  threshold  decreases  as  the 
number  of  data  snapshots  increases.  The  resolution  performance  of  MUSIC  as  a  function  of 
computational  load  is  compared  with  that  of  the  fast  projection  techniques,  DVSO  and  DVSO- 
COVAR,  in  Figure  23  and  26.  Figure  23  shows  A9^  versus  complex  multiplications  for  both 
MUSIC  and  DVSO-COVAR  techniques  for  an  array  of  10  elements.  The  DVSO  method 
performs  best  for  A9^  >  0, 12  BW.  The  DVSO-COVAR  technique  performs  best  for  0.09  BW 
<  A9^  <  0.12  BW.  Below  0.09  BW  the  MUSIC  technique  performs  best.  Figure  24  shows 
A9^  versus  complex  multiplications  for  an  array  of  40  elements.  For  A9^  >  0.13  BW,  the 
DVSO  technique  performs  best  while  for  A9aia  <  0.13  BW,  the  DVSO-COVAR  perform  best. 
As  K  increases,  the  region  in  which  DVSO-COVAR  performs  better  than  MUSIC  extends  to 
much  higher  resolution,  down  to  about  0.04  BW  for  a  40-element  array;  the  performance 
advantage  of  DVSO  over  MUSIC  in  the  intermediate  resolution  region  also  becomes  more 
pronounced.  It  should  be  noted  that  very  high  resolutions  (e.g.  0.04  BW)  are  usually  considered 
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Fig.  22  versus  S  for  MUSIC  with  SNR=40 
dB.  Monte-Carlo  simulations  with  100  trials. 


Complex  multipBcations 

Fig.  23  A0„u,  versus  complex  multiplications  for 
MUSIC  and  DVSO-COVAR  with  K=10  and 
SNR=40  dB.  Monte-Carlo  simulations  with  100 
trials. 
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Complex  multiplications 


Fig.  24  A0^  versus  complex  multiplications  for 
MUSIC  and  DVSO-COVAR  with  K=40  and 
SNR=40  dB.  Monte-Carlo  simulations  with  100 
trials. 
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to  be  not  practiced  for  MUSIC  because  of  sensitivity  to  calibration  errors  in  real  antenna 
systems.  The  sensitivity  of  DVSO  and  DVSO-COVAR  to  calibration  errors  has  not  been 
evaluated  and  is  beyond  the  scope  of  this  study. 


10.0  CONCLUSIONS 

We  have  compared  the  performance  of  three  projection  methods,  DVO,  DVSO  and 
DVSO-COVAR  with  the  performance  of  MUSIC  technique. 

The  selection  of  the  best  technique  for  a  particular  application  is  a  function  of  the 
computational  capability  available  and  the  degree  of  high  resolution  sought.  The  performance 
curves  presented  in  this  report  allow  a  more  nearly  optimum  trade-off  between  performance  and 
computational  load  than  has  heretofore  been  available.  Specifically,  if  a  moderate  degree  of 
super-resolution  is  required,  DVO  provides  the  best  combination  of  speed  versus  performance; 
a  somewhat  high  degree  of  resolution  leads  to  selection  of  DVSO;  still  higher  re'-''ution  requires 
DVSO-COVAR  with  its  high  computational  load.  The  above-mentioned  fast  projection 
techniques  provide  performance  superior  to  that  achievable  with  MUSIC  in  terms  of  resolution 
versus  computational  load  -  as  long  as  the  number  of  array  elements,  K,  is  high  compared  to 
the  number  of  jammers  L. 
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